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Abstract 

We investigate similarities in the micro-structural dynamics between externally driven and ac¬ 
tively driven nematics. Walls, lines of strong deformations in the director field, and topological 
defects are characteristic features of an active nematic. Similar structures form in driven passive 
nematics when there are inhomogeneities in imposed velocity gradients due to non-linear flow fields 
or geometrical constraints. Specifically, pressure driven flow of a tumbling passive nematic in an 
expanding-contracting channel produces walls and defects similar to those seen in active nematics. 
We also study the response of active nematics to external driving, confirming that imposed shear 
suppresses the hydrodynamic instabilities. We show that shear fields can lead to wall alignments 
and the localisation of active turbulence. 
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FIG. 1. Snapshots of the active turbulent state, (a) Flow field. The continuous lines are streamlines 
and the colour shading represents the vorticity field with blue and red shades respectively indicating 
clockwise and anticlockwise vorticity. (b) Corresponding nematic director held. Topological defects 
of charge ±1/2 are shown by red and blue markings respectively. Lines of strong deformations, 
which we shall term walls, are also distinctly visible in the director held. 


I. INTRODUCTION 

Active matter is composed of entities that produce their own energy and hence active 
systems are continually driven out of thermodynamic equilibrium [1-3]. Biological examples, 
across a range of length scales, include actin or microtubule hlaments powered by molecular 
motors [4, 5], bacterial suspensions [6, 7] and flocks of birds [8]. Moreover self-propelled 
particles can be fabricated to move by exploiting phoretic forces [9, 10]. Vertically vibrated 
granular rods are also often categorised as active systems [11]. 

Many active systems have nematic symmetry and can be described by the same hydrody¬ 
namic equations of motion as passive nematic liquid crystals, but with an additional active 
term in the stress tensor. The active stress generates a flow held whenever there is a gra¬ 
dient in the magnitude or direction of the nematic order. As a result, active nematics are 
hydro dynamically unstable, and dense active suspensions exhibit chaotic turbulent-like how 
states known as active turbulence [6, 7, 12]. Active turbulence is characterised by regions 
of strong vorticity in the velocity held as shown in Fig. 1(a). 
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Fig. 1(b) shows a typical corresponding director field. The hydrodynamic instabilities 
generate walls, lines of strong distortions in the director field [13, 14]. In an extensile active 
nematic, these walls are bend distortions as indicated in Fig. 1(b). These walls separate 
nematically ordered domains, and are sources of hydrodynamic stress that contribute to 
the chaotic flow fields shown in Fig. 1. However the walls are not stable structures. They 
disintegrate into one or more pairs of oppositely charged (±|) topological defects. Elastic 
forces and flow drive the defects apart, and they subsequently move around until oppositely 
charged defects meet and annihilate. Both defect creation and annihilation events regenerate 
nematic regions which again undergo instabilities leading to wall formation and providing a 
mechanism that sustains active turbulence. However, despite this qualitative understanding, 
a comprehensive theory describing the dynamics of walls and defects is yet to emerge. 

Active liquid crystals operate out of equilibrium because of their intrinsic activity which 
leads to an input of energy at every point in the fluid. However, nematics, both passive and 
active, can also be driven out of equilibrium by an applied shear or pressure gradient. In 
this paper we numerically solve the equations of motion of both active and passive nematics 
under shear and pressure-driven flow and compare the effects of external and internal driv¬ 
ing. In particular we show that, in certain circumstances, the walls and topological defects 
characteristic of active turbulence can also be a feature of driven, passive nematics. 

Firstly, in Secs. H and HI, we summarise the relevant hydrodynamic equations and the 
method that we use to solve them. In Sec. IV we discuss the effect of shear, confinement 
and Poiseuille flow on 2D passive nematics, concentrating in particular on the tumbling 
regime. We show that pressure driven flow of a tumbling nematic suspension leads to walls, 
which generate topological defects in a channel of non-uniform cross section. In Sec. V we 
turn to active nematics. Linear stability analysis has shown that there is a threshold shear 
beyond which active hydrodynamic instabilities are suppressed for aligning suspensions [15]. 
Solving the nonlinear equations we verify this analysis and show that this suppression of 
activity is observed in the tumbling regime as well. A corollary is that in Poiseuille flow, 
active turbulence is increasingly confined to the centre of the channel as the flow velocity 
increases. Sec. VI concludes the paper with a discussion of the relevance of our analysis of 
passively driven systems to active systems. 
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II. EQUATIONS OF MOTION 


The equations of nemato-hydrodynamics are described in this section. These are ap¬ 
plicable to both passive and active nematics, except for having an additional stress term 
when used to describe the latter. Therefore the system under consideration may be pas¬ 
sive, a molecular or colloidal liquid crystal, or active, such as a suspension of bacteria or 
microtubule bundles driven by molecular motors [1]. 

Both density and momentum are conserved variables, while orientational order is not. As 
well known for nematic liquid crystals, we use a second order traceless symmetric tensor Q 
to describe the orientational order. It is dehned as Q = |(3nn - I) where q is the largest 
eigenvalue of Q and is a measure of nematic degree of order, n is the director held and I 
is the identity tensor. We shall solve the equations of motion describing the evolution of 
momentum pu and order parameter Q [16, 17] 


p{Ot ~\~ 

(^dt ~\~ Uk^k^Qij tittup 


( 1 ) 

( 2 ) 


along with the equation of continuity. 

In addition to the advection, the response of Q to velocity gradients of the how held is 
also taken into account by the generalised advection term 


Sij (AT/j^ -|- ^ik') i^Qkj T ^kj/^') T i^Qik T ^ik/^')(^^Ekj ^kj') 

— 2\{Qij + 6ij/3){QkidkUi), (3) 


where Eij = {diUj + djUi)/2 is the strain rate tensor and Qij = {djUi — diUj)/2 is the vorticity 
tensor. 

Vorticity cj = V x u, the curl of velocity held describing the local rotation of the huid 
elements, is related to the vorticity tensor by eijkOJk = —2Qij where e is the Levi-Civita 
symbol. We will hnd that gradients of vorticity plays an important role in generating 
nematic director textures. 

The competition between elongation and rotation is determined by the alignment param¬ 
eter A which describes the aligning/tumbling motion of the director held in response to a 
simple shear how in both passive and active nematics. From a microscopic perspective A is 
determined by the degree of orientational order and the shape of the constituent particles. 


4 


The nematic suspension is usefully characterised in terms of Ai = (3g + 4)A/9g. For a tum¬ 
bling material, Ai < 1, the director undergoes continuous rotation in response to a simple 
shear flow. Instead, for Ai > 1 the director aligns at a hxed angle, known as the Leslie angle 
9i = A cos“^(l/Ai), with respect to flow direction [16-18]. 

We use the Landau-de Gennes free energy containing a gradient energy term and bulk 
energy terms, 

J- = K (dkQij)"^/2 AQijQji/2 C(QijQjiY/ 4, (4) 

to prescribe the equilibrium of the order parameter Q. Here, A, B and C are material con¬ 
stants that determine the conditions for the isotropic-nematic phase transition. We use the 
single elastic constant approximation where the elastic constant K assigns the same energy 
cost to splay, bend and twist deformations. The variational derivative of B with respect to 
Q defines the molecular field = —5B/5Qij -|- {5ij/?>)Ai{5B/5Qki)- The symmetric, trace¬ 
less matrix H describes the relaxation of Q in Eq. (2) at a rate controlled by the rotational 
diffusivity, D. 

The stresses acting on a fluid element may be classified into 

i) the viscous stress, = 2/iE'jj, 

ii) the elastic stress, 

= -PSij + 2A(g,, + 6,,/3){QkiHik) - XHMkj + 4,/3) 

Xi^Qik T ^ik/^')Hkj diQkl jrr\ ^ T QikHkj HikQkj j (h) 

OOjLJlk 

iii) the active stress, = —(Qij, 

where the modihed bulk pressure is P = pT — ^{dkQijY and the Newtonian viscosity of 
the suspension is p. The elastic stress term, generates ‘back-flow’, referring to flow 

helds generated in response to elastic stresses. 

The active stress has a strength ( which is proportional to the density of dipolar sources 
[19], unique to active systems. In passive systems C = 0. The sign of ( represents the 
symmetry of the dipolar flow held generated by the active units. In this paper, we only 
consider C > 0, corresponding to extensile systems, C < 0 corresponds to contractile systems. 
More details about these equations of motion and their application to passive and active 
systems can be found in [16, 17, 20-23]. 
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FIG. 2. Illustrations of various flows and boundary conditions studied in this work. The coordinate 
system used is also shown. The slanted lines indicate solid boundaries, (a) Lees-Edwards boundary 
conditions where the domain in the middle is sheared by its own images, (b) Nematic material 
conhned between two parallel plates with the top plate moving with a prescribed velocity, (c) A 
cubic velocity profile, u = is imposed, (d) A pressure gradient AP is imposed on the nematic 
material confined between two plates, (e) The channel in panel (d) is made non-uniform with an 
expanding and contracting cross section. 

The introduction of active stress to the nemato-hydrodynamic equations has non-trivial 
consequences. The nematic state is unstable to fluctuations of the order parameter field 
[19]. This results in spontaneous flows, which can be stabilised in a 1-D channel [22, 24]. In 
two dimensions the instability leads to active turbulence, shown in Fig. 1. 


III. NUMERICAL SOLUTION AND BOUNDARY CONDITIONS 

We use a hybrid lattice Boltzmann method on a D3Q15 lattice to numerically solve 
Eqs. (1) and (2) [22, 23, 25]. Density variations are considered to be negligible and the 
numerical scheme essentially yields the solution of the governing equations in the incom¬ 
pressible limit. We study the response of nematic liquid crystals to imposed driving by 
applying various boundary conditions of increasing complexity: 

i) Lees-Edwards boundary conditions: Typically a material is subjected to shear flow by 
confining it between two solid plates with one of the plates moving at a prescribed ve¬ 
locity. For nematic liquid crystals, anchoring on the solid boundary may effect the flow 
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and director profiles and therefore we implement Lees-Edwards boundary conditions, 
well known in molecular dynamics simulations, to apply a shear flow with no walls (see 
Fig. 2(a)). We shear the simulation domain by its own image, but with a larger mo¬ 
mentum. The director held, Q undergoes periodic boundary conditions perpendicular 
to the shearing direction, but with a modihed translational velocity. Details of the 
application of Lees-Edwards boundary conditions to lattice Boltzmann simulations are 
given in [26]. 

ii) Couette flow: We also subject the nematic material to a shear when conhned between 
two solid plates, with the top plate moving with a prescribed velocity and the bot¬ 
tom plate being stationary (see Fig. 2(b)). By specifying the director orientation on 
the plates, we can take into account the interaction of nematic material with solid 
boundaries. We choose homeotropic boundary conditions, with the director held on 
the boundaries constrained to lie perpendicular to the plates. Comparing the results 
to those using Lees-Edwards boundary conditions will help us to isolate the ehect of 
boundaries from that of shear. It turns out that this distinction is important in under¬ 
standing the wall formation. 

iii) Nonlinear velocity profile: We construct a model system to study the response of the 
director held to an imposed nonlinear how held. A method analogous to Lees-Edwards 
boundary conditions which generates a nonlinear velocity prohle is not available. There¬ 
fore, we subject the nematic material to a prescribed unidirectional how held with a 
cubic velocity prohle u = where the y direction is normal to the how (see Fig. 2(c)). 
We solve only Eq. (2) to give the response of the director to this how prohle, neglect¬ 
ing the mass and momentum conservation equations. This model system will help 
us to distinguish the ehect of nonlinear velocity prohles, which are also important in 
understanding the wall formation process without additional ehects arising from solid 
boundaries. 

iv) Poiseuille flow: A Poiseuille how prohle is produced by conhning the nematic mate¬ 
rial between two solid plates and subjecting it to an imposed pressure gradient (see 
Fig. 2(d)). As usual in lattice Boltzmann simulations, we impose a body force on the 
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fluid to generate the equivalent effect of a pressure gradient [27, 28]. No-slip velocity 
boundary conditions are applied on the bounding plates and homeotropic boundary 
conditions are specified for the director held. The resulting steady state velocity pro- 
hle is nonlinear. Thus we have created a natural situation to study the response of 
the nematic material to nonlinear velocity helds in the presence of solid boundaries, a 
combination of (ii) and (iii) above. 

v) Pressure driven flow in channel of non-uniform cross section: We consider Poiseuille 
how in a channel of nonuniform cross section along the imposed how direction (see 
Fig. 2(e)). Thus we introduce perturbations to the unidirectional how held. We consider 
a channel which is a periodic arrangement of expanding and contracting sections. The 
half period of the channel is chosen to be 100 lattice points and the expanding and 
contracting angles are 10°. As before no-slip/ homeotropic boundary conditions are 
imposed on the plates for the velocity/director helds. 

The parameters used in the simulations are D = 0.34, A = 0.0, B = —0.3, C = 0.3, 
K = 0.01, C = 0.01, /i = 2/3, in lattice units, unless mentioned otherwise. Simulations 
were performed in a domain of size 100 x 100 except in case (v) above where a domain 
of 200 X 100 was chosen. The choice of boundary conditions for the director held on the 
solid wall, namely homeotropic or planar anchoring, does not ahect the conclusions of this 
work and therefore we chose to report results of simulations using the former alone. The 
imposed shear rate (with or without solid boundaries) is 10“^. The applied pressure gradient 
generating Poiseuille how in the channel geometries is 5 x 10“®. As usual in lattice Boltzmann 
schemes, discrete space and time steps are chosen as unity and therefore all quantities can 
be converted to physical units in a way which depends on the material of interest [29-31]. 
This procedure will demand choosing appropriate mass, length and time scales. Molecular 
liquid crystals are typically ~ lOnm in length, exhibit a rotational dihusivity ~ 100-1000 
(Pa-s)“^ and have elastic constants ~ IpN [16]; and these can be used as an appropriate 
basis for conversion to physical units. However, for active liquid crystals, similar data are 
not available yet. A similar system may be FD viruses, passive colloidal liquid crystals, 
which are microns in length, have rotational dihusivity ~ 20 (Pa-s)“^ and elastic constants 
IpN. 
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FIG. 3. Response of the nematic director field to a simple shear flow imposed by Lees-Edwards 
boundary conditions (see Fig. 2(a)). The director held (only 9 grid points are shown here), which 
rotates in the clockwise direction following the vorticity, is shown at equal intervals of time, colour- 
coded blue to magenta (see text for a detailed explanation), (a) Ai > 1: the director held aligns 
at the Leslie angle, (b) Ai < 1: the director tumbles. The angular velocity is largest when the 
director is along the how gradient direction. 


IV. DRIVEN PASSIVE NEMATICS 

We first summarise the response of a passive, two-dimensional, nematic suspension to a 
simple shear flow imposed using Lees-Edwards boundary conditions. Fig. 3 compares the 
behaviour of aligning and tumbling materials. For Ai > 1 the director aligns at the Leslie 
angle, 6i = ^ cos“^(l/Ai), with respect to the flow direction as shown in Fig. 3(a). For 
Ai < 1 the director undergoes continuous tumbling, and Fig. 3(b) shows its time evolution 
during part of a rotation period. This, and subsequent, hgures should be read as follows: 
The time has been given a color coding represented by the colorbar on the right hand side. 
The orientation of the director at the initial time is shown in blue. As time proceeds, the 
director rotates following the imposed vorticity. This rotation is in the clockwise direction 
as indicated by the curved arrow. The hnal orientation of the director is shown in magenta. 
Intermediate orientations at equally spaced time intervals, which correspond to the color- 
bar, are also shown. Just as for a single particle in a simple shear flow, the maximum angular 
velocity occurs when the director held is oriented along the gradient direction (y-axis) and 
the minimum when it is in the direction of the how (x-axis). We also give the diherence in 
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FIG. 4. Wall formation in the director field when an (a) flow aligning and (b) flow tumbling nematic 
is subjected to Couette flow. Only the top half of the domains are shown. The illustrations on 
the LHS of panels (a) and (b) follow the conventions of Fig. 3. The RH plot is the director field 
corresponding to the last frame of the LH plot. Imposed homeotropic boundary conditions on the 
plates do not allow spatially homogeneous tumbling of the director field. Thus, in the flow aligning 
case (a), a half wall is formed near the top solid plate and the director aligns at the Leslie angle 
everywhere else. In (b) the director field tumbles, and spatial inhomogeneities generate walls. The 
flow profiles in each case are also shown. 


angle between the chosen final and initial time - this will be relevant when the number of 
rotations varies with position. 


So far we have observed that nematics subjected to a simple shear without any solid 
boundaries and confined to two dimensions remain homogeneous in space. There are no 
walls or topological defects that are characteristic features of active turbulence. We now 
investigate the essential ingredients that are required for walls and topological defects to 
form in passive nematics. 


10 









752 -^ 

755 -^ 

745 ^ 

745 --' —^ 

735—" —"" 

68 >^ ^ 
62 ^^ ^ 
54V> ^ 

45° ^ ^ 




(a) 


(b) 


FIG. 5. Wall formation in the director field when an (a) flow aligning and (b) flow tumbling 
nematic is subjected to a cubic velocity profile. The illustrations on the LHS of panels (a) and (b) 
follow the conventions of Fig. 3. The RH plot is the director field corresponding to the last frame 
of the LH plot. The spatial gradient in vorticity causes the director field to rotate with different 
angular velocities giving (a) an inhomogeneous orientation field in the aligning case and (b) walls 
in a tumbling material. The flow profiles in each case are also shown. 


A. Walls in the director field 


Both walls and topological defects are inhomogeneous director helds. Therefore we now 
investigate whether inhomogeneous forcing can generate wall-like structures. Indeed we hnd 
that walls naturally form if the angular velocity of the director held in response to how 
gradients or elastic forces is inhomogeneous in space. We show that there are two diherent 
ways to impose director inhomogeneities (i) via the introduction of solid boundaries and (ii) 
by imposing nonlinear how helds. 

We hrst consider Couette how with homeotropic anchoring specihed on both bounding 
plates. The simulations were started with director held normal to the plates, consistent with 
the boundary conditions. The top plate is allowed to move and the resulting conhgurations 
of the director held are illustrated in Fig. 4 for both a shear aligning and a shear tumbling 
material. Only the top half of the channel is shown, as the conhgurations have rehection 
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symmetry about the centre of the channel. For flow aligning nematics, Fig. 4(a) shows 
that the bulk of the sample is rotated by approximately 74° (in this case) and aligned 
with the flow direction at the Leslie angle (~ 16°). This is the same behaviour as for 
Lees-Edwards boundary conditions. However, near the solid boundary, the homeotropic 
anchoring competes with the flow aligning tendency of the director held, thus generating 
a smooth variation in the orientation. The result is the formation of a half wall near each 
boundary. 

In the case of tumbling nematics in Couette how there is a similar competition between the 
anchoring boundary conditions and the rotation of the director held as shown in Fig. 4(b). 
Unlike the aligning suspensions, the response of tumbling materials is not uniform through¬ 
out the bulk. This is because, as time proceeds, the director held continues to rotate instead 
of approaching a constant angle of orientation and the total angle rotated increases between 
the edges and the centre of the channel (from approximately 21° to 429° in this case). This 
spatially diherential tumbling results in walls in the director held as shown in Fig. 4(b). We 
stress that Fig. 4(b) is not a steady state, tumbling of the director held continues in response 
to the imposed how and how gradients. Thus a dynamic situation persists with continual 
formation and destruction of walls in the director held. 

Secondly, we show that diherential rotation rates can also be generated by imposing 
nonlinear velocity helds. We solve Eq. (2) with an imposed cubic velocity prohle. As before 
the simulations were started with an initial director held oriented along the y direction. The 
results are illustrated in Fig. 5. 

In the case of how aligning nematics, the director held rotates towards the Leslie angle. 
Note that the Leslie angle is independent of the shear rate and depends only on the value of 
Ai. The rotation is fastest where the velocity gradient is largest (near top of the domain). 
The bottom of the domain experiences smaller shear and had not rotated enough to reach 
the Leslie angle during this simulation. In the case of tumbling nematics, the director held 
continuously rotates. Due to the cubic velocity prohle, vorticity is not constant in space 
and hence diherent regions along the y direction rotate with diherent angular velocities. 
While the director at the top of the domain has rotated more than a period, the bottom has 
rotated only by a quarter of a period. The result is the formation of walls. Again, Fig. 5(b) 
is not a steady state picture, but a representation of a temporally evolving state. 

Thus the diherential rotation of the director held in response to vorticity gradients, caused 
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FIG. 6. Wall formation in the Poiseuille flow of nematic liquid crystals. The illustrations on the 
LHS of panels (a) and (b) follow the conventions of Fig. 3. The RH plot is the director field 
corresponding to the last frame of the the LH plot, (a) Aligning suspensions: the solid boundaries 
generate half walls and flow gradients generate a wall at the centre of the channel, (b) Tumbling 
suspensions: the solid boundaries and nonlinear velocity profile act together to generate walls. The 
flow profiles in each case are also shown. 


by boundaries or varying shear, generates inhomogeneities in the director field which, in tum¬ 
bling suspensions, may be identified as walls. In a Poiseuille flow geometry both boundaries 
and shear gradients are present. Therefore we now consider the flow of a nematic material 
between two solid plates, with homeotropic boundary conditions on the director, under a 
pressure gradient. The results are shown in Fig. 6. The simulations were started with an 
initial director field in the vertical direction. 

For aligning nematics, the director field rotates and reaches a steady state at the Leslie 
angle everywhere except near the boundaries and in the centre of the channel. As for Couette 
flow, the competition between the flow alignment and the boundary anchoring leads to the 
formation of half-wall structures near both boundaries. In the centre of the channel, where 
the slope of the velocity profile changes sign, the director held rotates to generate a wall. As 
the driving pressure gradient increases, the bend conhguration of this wall (known as a ‘bow’ 
mode in the literature) is replaced by a ‘peak’ conhguration [32]. Such changes in director 


13 






















configuration are known to change the rheology of the material [33]. Flow of an aligning 
nematic under an imposed pressure gradient with various anchoring conditions has been 
recently investigated in [34] due to its importance in microfluidics. This study also reports 
similar director conhgurations. In the case of pressure driven flow in three dimensional 
channels, the director held need not be constrained to two dimensional planes and walls can 
be replaced by three dimensional director conhgurations arising from a competition between 
how and anchoring ehects [35]. Small diherences in anchoring are known to lead to large 
diherences in the director conhgurations [36]. 

In case of how tumbling nematics under pressure driven how we observe several walls gen¬ 
erated by the boundary anchoring and the nonlinear velocity prohle acting together. Again, 
the director held does not reach a steady state, but instead walls are continually created, 
destroyed and displaced. The corresponding velocity prohles, shown for both aligning and 
tumbling nematics in Fig. 6, have no noticeable structure at the location of walls. 

The complex behaviour arising from the competition between how and ordering ehects 
in passive nematic liquid crystals and colloidal rods have been experimentally demonstrated 
in some recent investigations [37-40]. Our simulations are distinct from these experiments 
and most previous work in that we constrain the director held to two dimensions and do not 
allow the director to escape in the vorticity direction on the application of a shear force. This 
is to give a direct comparison with experimental active systems such as actin or microtubule 
hlaments driven by molecular motors [5] and bacterial suspensions [41] where the active 
nematic material was constrained in two dimensional hlms. Moreover, in these cases, the 
length characterising the constituents is on much larger scales than typical molecular liquid 
crystals making escape into the third dimension less likely. 


B. Defect formation from walls 

So far, we have been looking at the formation of walls in nematics. For all the diherent 
situations that we have studied, walls did not show any instabilities leading to defect for¬ 
mation as observed in active turbulence (see Fig. 1(b)). This is a consequence of the one 
dimensional nature of the imposed forces. Walls can become unstable if there are sufficiently 
large perturbations or inhomogeneities. We introduce the latter by using a channel of non- 
uniform cross section as shown in Fig. 7. We consider tumbling nematics in Poiseuille flow 
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(a) (b) 

FIG. 7. (a) Tumbling passive nematics subject to pressure-driven flow in a channel of non-uniform 
cross section generate walls and defects, (b) corresponding flow profile with the average velocity 
in X—direction subtracted. Both streamlines and vorticity field (colour shading) are shown. 


with homeotropic anchoring at the solid boundaries. As discussed in the previous section 
walls are generated due to shear. However, now the walls are slightly perturbed along the 
flow direction due to variations in the flow and director held in the expanding-contracting 
cross section of the channel. Such perturbations are sufficient to drive the walls unstable 
and they break down to form a pair of oppositely charged topological defects. The defects 
both destroy the walls and regenerate nematic orientational order. Due to the imposed how, 
they are carried away and execute their own motion driven by elastic forces and back how 
ehects, annihilating on hnding defects of opposite charge. Thus, walls and defects are con¬ 
tinuously created and destroyed in a tumbling nematic suspension when driven in a channel 
of non-uniform cross section. 

This can be considered the counterpart of active nematics where these phenomena also 
occur but driven by active forces. However, unlike for active turbulence the presence of walls 
and defects in passive tumbling nematics does not alter the how held substantially. Panel 
(b) of Fig. 7 shows the vorticity and streamlines of the how held once the average horizontal 
velocity is subtracted. The circulating hows are a consequence of mass conservation in the 
expanding-contracting cross section of the channel and are not related to walls or defects. 
We also note that walls can disintegrate into topological defects even in a channel of uniform 
cross section if thermal huctuations are sufficiently large. 
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V. DRIVEN ACTIVE NEMATICS 


A. Suppression of hydrodynamics instabilities under shear 

We now consider driven active nematics, thus studying situations in which both active 
and external forces are acting together. Firstly, active nematics are subjected to simple 
shear flow by applying Lees-Edwards boundary conditions. This situation was studied by 
Muhuri et al. [15] who used linear stability analysis to show that hydrodynamic instabilities 
are suppressed by shear and that the critical shear rate needed to stabilise the active nematic 
increases with decreasing aligning parameter, A. They obtained a phase diagram of Pe^ vs 
A where Pe^ = 2^'y/C, measures the ratio of the imposed driving to the active driving, with 
7 the imposed shear rate. We construct a similar plot. Fig. 8(a), by numerically solving 
the entire set of nonlinear governing equations, conhrming their analysis beyond the linear 
regime. Moreover our simulations show that as A decreases into the tumbling regime the 
critical shear rate for stability continues to increase, consistent with the fact that tumbling 
materials have an intrinsic tendency to generate walls. However, simulations show that when 
A ~ 0.5 the critical shear rate decreases again. The reasons for this change in behaviour are 
not yet clear. 

Active nematics exhibit an interesting behaviour at shear rates close to the value required 
for stabilising the nematic. The walls generated by the active stresses themselves respond 
to the imposed shear field. Behaving effectively as a single structure, they align close to the 
flow direction. This wall alignment regime is illustrated in Fig. 8(b). Defects may continue 
to be produced, so that the walls continually disappear, and then re-form. 

B. Poiseuille flow 

The competition between the stabilising effect of shear and the destabilising effect of 
the activity is also seen in an active nematic under Poiseuille flow. The imposed pressure 
gradient generates a velocity prohle with a maximum at the centre of the channel as in 
Fig. 6. Shear rates are smallest at the centre of the channel and thus active stress can 
dominate the imposed shear rate and walls and defects form. Nearer the walls the shear 
rates are sufficiently large that active hydrodynamic instabilities are suppressed and the 
nematics align at the Leslie angle. Hence active turbulence is confined to the centre of the 
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(a)C = 0.002 


(b)C = 0.002, A = 0.4 



FIG. 8. Suppression of the hydrodynamic instabilities of an active nematic by shear, (a) phase 
diagram showing the stable and unstable regions as a function of Pea, the ratio of imposed driving 
to active driving, and A, the alignment parameter. Simulations were performed in a domain of size 
100 X 200. Both aligning and tumbling regimes are shown, (b) alignment of walls due to simple 
shear before the hydrodynamic instabilities are completely suppressed, (c) active turbulence in 
Poiseuille flow. Instabilities and thus active turbulence are seen only at the centre of the channel 
where the shear rate due to the imposed pressure gradient is weak, (d) The width of the turbulent 
layer increases as activity increases. 
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(b) 

FIG. 9. Tumbling active nematics driven by an imposed pressure gradient in a non-uniform channel, 
(a) walls and defects, (b) flow field with the average velocity in the x—direction subtracted. Both 
streamlines and vorticity (colour-shading) are shown, (c) Number of defects as a function of activity 
(. The number of defects tends to a constant value as the activity is reduced. 

channel as illustrated in Fig. 8(c). If activity is further increased the width of the turbulent 
layer in the channel is increased. At sufficiently large activities, activity driven flow will 
dominate the imposed flow. 

To complete the picture and to compare to passive nematics we now consider the be¬ 
haviour of active tumbling nematics in a periodically expanding-contracting channel driven 
by a pressure gradient. The results for active nematics are illustrated in Fig. 9 (compare 
Fig. 7 for the passive case). Panel (a) shows the director field, dominated by walls and 
defects as both the active stress and the imposed pressure gradient generate these director 
textures. Panel (b) shows the corresponding vorticity and stream lines, again after subtract¬ 
ing the mean flow along the channel. This is very different to the passive case as the defects 
produce additional flows through the active term in the stress tensor. 

To link the active and passive systems we plot the number of defects at different activity 


18 




















(a) (b) 

FIG. 10. (a) Spatial correlations of the nematic director field for various values of the alignment 
parameter A for an active nematic (( = 0.01). Both aligning and tumbling regimes are shown. 
The characteristic length of the nematic order decreases as A decreases due to the appearance of 
more walls and topological defects. This is also apparent from (b) the number of defects as a 

function of the alignment parameter A. 


levels in Fig. 9(c). As activity decreases, so the number of defects reduces. In the case of a 
non-driven system, the number of defects would reduce to zero at sufficiently small activities 
[31]. However, here, a hnite number of defects is generated even at zero activity due to the 
driving alone. 


VI. DISCUSSION 

Our simulations were in 2D with the director held conhned to the plane to give a direct 
comparison to most previous work on active systems. Therefore we see how aligning and 
tumbling regimes under shear but not the instabilities that lead to kayaking, log rolling or 
rheochaos which require escape of the director into the third dimension, or wagging which 
would only appear at much higher shear rates [18, 42-44]. Indeed, little is known about 
active turbulence in 3D, even in the absence of shear. 

By analysing the response of passive nematics to external driving forces, we have found 
that spatial inhomogeneities in the angular velocity of the director held, either due to elastic 
forces blocking the rotation or spatially dihering vorticity helds, gives rise to walls. Thus 
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tumbling nematics have an intrinsic tendency to generate more walls than aligning nematics. 

Our finding that vorticity gradients help in the formation of walls is relevant to active 
turbulence. In active turbulence the walls and topological defects that are formed are the 
sources of vorticity, which then diffuses generating vorticity gradients. This helps in the 
formation of new walls which then sustain the active turbulence. To demonstrate this 
we analysed active turbulent patterns as a function of the alignment parameter A. Our 
observations are quantified in Fig. 10. In panel (a) the normalised spatial correlations of the 
nematic director field (Q(0) ; Q{R)) / (Q(0) : Q(0)) are plotted as a function of the radial 
distance R for a range of values of A. As A decreases, there is a clear reduction in the 
characteristic correlation length of the the director field. This indicates that more walls are 
formed as the tumbling tendency of the director field increases. The increased number of 
walls gives rise to more topological defects, as apparent in Fig. 10(b). 

To conclude, we have shown that in passive nematics walls can be generated due to bound¬ 
aries or imposed nonlinear flow fields. The effects are more pronounced when the nematic 
is tumbling. The walls can be easily destabilised by perturbing them, say with channels 
of non-uniform cross section. Then walls disintegrate into topological defects producing a 
director field analogous to that observed in active turbulence. Subjecting active nematics to 
external driving invokes a direct competition between active forces and externally imposed 
forces. When the imposed driving is large, active nematics behave like passive nematics. 
When the forces are comparable the competition leads to wall alignment in a simple shear 
flow and to active turbulence being confined to the centre of the channel in Poiseuille flow. 
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